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Abstract 

In this paper, we calculate the branching ratios and CP-violating asymmetries for — > 
fo(980)K, /o(1500)-KT within Perturbative QCD approach based on kx factorization. If the mix- 
ing angle 9 falls into the range of 25° < 9 < 40°, the branching ratio of B® — > fo(980)K is 
2.0 x 10~ 6 < B{B Q S -> /o(980)iO < 2.6 x 10" 6 , while 9 lies in the range of 140° < 9 < 165°, 
B(B° S -> /o(980)A") is about 6.5 x lO" 7 . As to the decay B{B° S -> / Q (1500)if), when the mixing 
scheme | /o(1500)) = 0.84 | ss) — 0.54 | nn) for /o(1500) is used, it is difficult to determine which 
scenario is more preferable than the other one from the branching ratios for these two scenarios, 
because they are both close to 1.0 x 10 -6 . But there exists large difference in the form factor 
i^^s-^/oCiSoo) £ or £ WQ scenar ios. 
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I. INTRODUCTION 



For scalars' mysterious structure, it arose much interest in both theory and experiment. 
In order to uncover the inner structures, many approaches are used to research the B u ^ 
decay modes with a scalar meson in the final states, such as the generalized factorization 
approach (l^QCD factorization approach (QCDF) Perturbative QCD (PQCD) 

approach |a-llO|| . But as to B® meson, these decay modes haven not been well studied by 
theory. The role of scalar particles in B° s decays should be given much more noticeable, 
because analyses of the corresponding decays can also provide a unique insight to the 
mysterious structure of the scalar mesons. Here we we will study the branching ratios 
and CP asymmetries of B® — > fo(980)K, /o(1500)ir within Perturbative QCD approach. 
The fundamental concept of this approach is factorization theorem, which states that 
the nonperturbative dynamics pross can be separated from a high-energy QCD process. 
The remaining part, being infrared finite, is calculable in perturbation theory. So a full 
amplitude is expressed as the convolution of perturbative hard kernels with hadron wave 
functions. There is a parton momentum fraction x in the former, both x and in the 
latter. Because x must be integrated over in the range between and 1, the end-point 
region with a small x is not avoidable. If there is a singularity developed in a formula, 

factorization should be employed [ll, Here k? denotes parton transverse mometa. 
A wave function, because of its nonperturbative origin, is not calculable, but process 
independent. So it can be determined by some means, such as QCD sum rules and lattice 
theory or extracted from experimental data. On the experimental side, some of B® decays 
involved a scalar in the final states might be observed in the Large Hadron Collider beauty 



experiments (LHC-b) [13|, [14J. In order to make precision studies of rare decays in the B- 
meson systems, the LHC-b detector is designed to exploit the large number of b-hadrons 
produced. Furthermore, it can reconstruct a B-decay vertex with very good resolution, 
which is essential for studying the rapidly oscillating B s mesons. So the studies of these 
decay modes of B® are necessary in the next a few years. 

It is organized as follows: In Sect JTT], we introduce the input parameters including 
the decay constants and light-cone distribution amplitudes. In Sec JIIIl we then apply 
PQCD approach to calculate analytically the branching ratios and CP asymmetries for 
our considered decays. The final part contains our numerical results and discussions. 



II. INPUT PARAMETERS 

For the underlying structure of the scalar mesons is still under controversy, there 
are two typical schemes for the classification to them (lol. [l6|. The scenario I (SI): the 
nonet mesons below 1 GeV, including f (600), fo (980), K* (800) and a (980), are usu- 
ally viewed as the lowest lying qq states, while the nonet ones near 1.5 GeV, including 
/o(1370),/o(1500)// (1700),A:*(1430) and a (1450), are suggested as the first excited 
states. In the scenario II (SB), the nonet mesons near 1.5 GeV are treated as qq ground 
states, while the nonet mesons below 1 GeV are exotic states beyond the quark model 
such as four-quark bound states. In order to make quantitative predictions, we identify 
/o(980) as a mixture of ss and nn = (uu + dd)/y/2, that is 

|/ (980)) = \ss)cos6+ \nn) sin6, (1) 
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where the mixing angle 9 is taken in the ranges of 25° < 9 < 40° and 140° < 9 < 165° 17 



Certainly, / (1500) can be treated as a qq state in both SI and SII. We considered that 
the meson / (1500) and / (980) have the same component structure but with different 
mixing angle. 

For the the neutral scalar mesons /o(980) and /q(1500) cannot be produced via the 
vector current, we have (fo(p)\q2l^qi\0) = 0. Taking the mixing into account, the scalar 
current (/oG°)| ( Z2<?i|0) = msfs can be written as: 

(mdd\0) = (muu\0) = ^rn f Jf , (f»\85\0) = m f J s fo , (2) 

where f^ 1 '^ represent for the light-cone distribution amplitudes for nn and ss components, 
respectively. Using the QCD sum rules method, one can find the scale-dependent scalar 
decay constants /? and /? are very close 0]. So we shall assume /J? = fj o and denote 
them as ff in the following. 

The twist- 2 and twist-3 light-cone distribution amplitudes (LCDAs) for different com- 
ponents of /o are defined by 

(fa(p)\q(z)iq(0) s \0) = J dxe^ z {^ h {x) + m fo $ S fo (x) + m f M + i- ~ 

(3) 

here we assume /^(p) and /o(p) are same and denote them as fo(p), n + and n_ are light- 
like vectors: n + = (1,0,0^),^- = (0,1, 0^). The normalization can be related to the 
decay constants 



1 dx* k (x) = / dx*l(x) = 0, / dx*%(x) = (4) 



The wave function for K meson is given as 



Sjc(P, x, C) = ^ip75 [^fttix) + m«$ p K (x) + Qm^iU - v ■ n)<S> T K {x)] , (5) 
v ^ iV c 

where P and x are the momentum and the momentum fraction of K meson, respectively. 
The parameter ( is either +1 or —1 depending on the assignment of the momentum 
fraction x. 

In general, the £L meson is treated as heavy-light system and its Lorentz structure can 



be written as 18, 19 



$ Bs = ^=tffl. + M fl .)7te^.(*i)- (6) 



For the contribution of <pB a is numerically small [20| and has been neglected. 



III. THEORETICAL FRAMEWORK AND PERTURB ATIVE CALCULA- 
TIONS 

Under the two-quark model for the scalar mesons supposition, we would like to use 
PQCD approach to study B° s — > f {9S0)K, f (1500)K decays. In this approach, the decay 
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amplitude is separated into soft, hard, and harder dynamics characterized by different 
energy scales (i,m_B s , M w ). It is conceptually written as the convolution, 



A(B° S -> f K) ~ Jd 4 hd 4 k 2 d 4 k 3 Tv[C(t)^ Bs (ki)^f (k 2 )^K(h)H(k 1 ,k 2 ,h,t)}, (7) 

where k^s are momenta of anti-quarks included in each meson, and Tr denotes the trace 
over Dirac and color indices. C (t) is the Wilson coefficient which results from the radiative 
corrections at a short distance. The function H(ki, k 2 , k 3 , t) describes the four quark 
operator and the spectator quark connected by a hard gluon whose q 2 is in the order of 
AMs a , and includes the 0(y AMr s ) hard dynamics. Therefore, this hard part H can be 
perturbatively calculated. 

Since the b quark is rather heavy, we consider the B° s meson at rest for simplicity. It 
is convenient to use light-cone coordinate {p + ,p~,Px) to describe the meson's momenta, 



V 



u ±p d ), and p T = (p 1 ,P^ 



V2 



(8) 



Using these coordinates the B° s meson and the two final state meson momenta can be 
written as 



Mr Mr 

* 1,1, Or , P 2 = ^(1,0,0 T ), P 3 



Mr 



y/2 



s/2 



s/2 



(0,1, T ), 



(9) 



respectively. The meson masses have been neglected. Putting the anti-quark momenta in 
B° s) /o and K mesons as k\, k 2 , and k 3 , respectively, we can choose 

k 1 = (x 1 Pf,0,V: 1 T), k 2 = (x 2 P+,0,k 2T ), k 3 = {0,x 3 P 3 ,k 3T ). (10) 

For our considered decay channels, the integration over fcf, k 2 , and k 3 in equation (JTj) 
will lead to 



A(B S foK) ~ J dx\dx 2 dx 3 b\db]b 2 db 2 b 3 db 3 

•Tr [C(t)$ Bs (x l7 6i)$ /0 (x 2 , b 2 )<S> K (x 3 , b 3 )H(x l} b h t)S t (x t ) e~ s(t) ] , (11) 

where 6j is the conjugate space coordinate of kiT, and t is the largest energy scale in 
function H(xi,bi,t). The large double logarithms (Id 2 Xi) on the longitudinal direction 
are summed by the threshold resummation 21], and they lead to S t (xi), which smears 
the end-point singularities on x«. The last term e~ s ^ is the Sudakov form factor which 
suppresses the soft dynamics effectively [22|. Thus it makes the perturbative calculation 
of the hard part H applicable at intermediate scale, i.e., M# s scale. 

We will calculate analytically the function H(xi,bi,t) for B° s — > f$K decays in the 
leading-order and give the convoluted amplitudes. For our considered decays, the related 
weak effective Hamiltonian H e ff can be written as [23[ 



H. 



eff 



71 



q=u,c 



10 



{CMO q M + c 2 (ti)o«M) ]T go*) ofr) 



i=3 



(12) 
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with the Fermi constant Gf = 1.16639 x 10 5 GeV 2 , and the CKM matrix elements V. 
We specify below the operators in % e fj for b — > d transition 

0\ = da'yVLup ■ u^f.Lba , O2 = d_ a ^Lu a • U^^Lbp , 

o 5 = d a YLb a ■ J2 q < q'pinRq'p , o 6 = d a ^Lb p ■ J2 q < q'pinRq'a » (13) 

7 = \d a ^Lb a ■ 2 g / e^q'^Rq^ , 8 = \d a ^Lbp ■ Y^ q , e^q'^^Rq^ , 
9 = \d a ^Lb a ■ Yjq> Zq'q'plvLq'p , O 10 = \d a ^Lb p ■ J2 q > Zq4f\ 'i* L( l' a > 

where a and (3 are the SU(3) color indices; L and R are the left- and right-handed 
projection operators with L = (1 — 75), R = (1 + 75). The sum over q' runs over the quark 
fields that are active at the scale \i = 0(m,B 3 ), i.e., (q'e{u, d, s, c, b}). 

There are eight type diagrams contributing to the B s — > f Q K decays are illustrated 
in figured) For the factorizable emission diagrams (a) and (b), Operators 01,2,3,4,9,10 are 
(V — A)(V — A) currents, and the operators 05,6,7,8 have a structure of (V — A)(V + A), 
the sum of the their amplitudes are written as F e f and F^}. In some other cases, we 
need to do Fierz transformation for the (V — A) (V + A) operators and get (S — P)(S + P) 
ones which hold right flavor and color structure for factorization work. The contribution 
from operator type (S — P)(S + P) is written as F^?\ Similarly, for the facorizable 
annihilation diagrams (g) and (h), The contributions from (V — A) (V — A), (V — A)(V + 
A),(S — P)(S + P) these three kinds of operators are F a f ,F^ and F^, respectively. 
For the nonfactorizable emission (annihilation) diagrams (c) and (d) ((e) and (f)), these 
three kinds of contributions can be written as M e ( a )j , M^K^, M^A^, respectively. Since 

these amplitudes are similar to those for the decays B — > / (980)iT(7r, 77^) @, Si or 
B — > 00(980)^ [To}, we just need to replace some corresponding wave functions and 
parameters. It is the same with the amplitudes for the fo and K° exchanging diagrams. 

Combining the contributions from different diagrams, the total decay amplitudes for 
these decays can be written as 

M{B° S f K) = M ss xcos9 + l=M n nsme, (14) 

a/2 



with 

Mss = -VM 



1 x , „DO „DO> , 1 

2* 



(F ef0 + F af0 )(a 4 - -a w ) + (F e f o + F^)(o 6 - ^s) 



1 - * rPl , n/rPl 1 



+(M e/0 + M afo )(C 3 - -C 9 ) + (Mff + M£ )(C 5 - -Cr) 



+M eK (C 4 - ~C W ) + M™(C 6 - l -C 8 



(15) 



M r 



V ub V: d M eK C 2 - VM 



(F% + F%)(a e - -a 8 ) + {M% + M%)(C S 



+M eK (C 3 + 2C 4 + i(C 10 - C 9 )) + M e ^(2C 6 + l -Cz) + M aK (C 3 - l -C 9 ) 
+F aK (a 4: — -aio) 



(16) 



5 




fo 



(a) 



Be 



K° 



(c) 



f m 



(b) 



K° 



(d) 




-0 




K 



fo 



A 





FIG. 1: Diagrams contributing to the B® — > foK° decays . 



where the combinations of the Wilson coefficients are defined as usual 24, 25 



a x = C 2 + d/3, a 3 = C 3 + C 4 /3, a 5 = C 5 + C 6 /3, a 7 = C 7 + C 8 /3, a 9 = C 9 + C 10 /3, 
0,2 = Ci + C2/3, ct4 = C4 + C3/3, Q-6 = C*6 + C5/3, Ct8 = C*8 + Cj/3, <2io = C\q + Cg/3. 



;i7) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

The twist-2 LCDA can be expanded in the Gegenbauer polynomials 

$ /o O,/i) 



^-i=/ /o ( M )6x(l - x) £ iU^ /2 (2a; - f), 



(19) 



m=l 



where B m (jj) and Cm 2 (a;) are the Gegenbauer moments and Gegenbauer polynomials, 
respectively. The values for Gegenbauer moments and the decay constants are taken (at 
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TABLE I: Input parameters used in the numerical calculation [31. \26\. 



Masses 


"Vo(980) = 0.98 GeV, 




= 1.7 GeV, 




"Vo(i500) = 1-5 GeV 


M Ba 


= 5.37 GeV, 


Decay constants 


f Bs = 0.23 GeV, 


Sk 


= 0.16 GeV, 


Lifetimes 


t b o = 1.466 x 10" 12 s, 


CKM 


V tb = 0.9997, 


V u = 


0.0082e- i2L6 °, 




V ud = 0.974, 


v ub = 


0.00367e- i60 ° . 



scale /j = lGeV) as @, Q 



Scenario I : // (g 8 o) 
B x (980) 
Si(1500) 
Scenario II : // (i 5 oo) 
Si(1500) 



(0.37 ± 0.02)GeV, // o( i 5 oo) = -(0.255 ± 0.03)GeV, 

-0.78 ±0.08, £ 3 (980) = 0.02 ±0.07, 

0.80 ±0.40, ^3(1500) = -1.32 ±0.14; 
(0.49 ± 0.05)GeV, 

-0.48 ±0.11, ^3(1500) = -0.37 ±0.20. (20) 



As for the explicit form of the Gegenbauer moments for the twist-3 distribution ampli 
def 
form 



tudes $j o and $ J o , they have not been studied in the literature, so we adopt the asymptotic 



2y/2N ( 



Jfo, $ 



2y/2N ( 



J f0 (l-2x). 



(21) 



In the each twist LCDAs, the appearance of Gegenbauer polynomials is from the 
expansion of non-local operator into local conformal operators 23]. Sometimes, the twist- 
3 contributions are important, because they can be enhanced due to some mechanism like 
chiral enhancement, especially in the condition of the leading twist contributions being 
small or zero. 

The twist-2 kaon distribution amplitude and the twist-3 ones <3>^ and <3>^ have 
been parametrized as 



K 



2y/2N c 
fK 



2y/2N c 

fK 



2V2N, 

The meson's wave function can be written as 



6x(l - x) [1 + 0.51(1 - 2x) + 0.3(5(1 - 2x) 2 - 1 

[1 + 0.24(3(1 - 2a;) 2 - 1) 

-0.12/8(3 - 30(1 - 2x) 2 + 35(1 - 2xf ' 
(1 - 2x) [I + 0.35(1 - lOx + 10x 2 )] . 



(22) 

(23) 
(24) 



M 2 x 2 1 

•> B .(x,b) = N Ba x 2 (l - X ) 2 expi-^r- - -(u bs b)% 



(25) 



where Ub B is a free parameter and we take u\ )a = 0.5 ±0.05 GeV in numerical calculations, 



and Nb s = 63.67 is the normalization factor for Ub a = 0.5 [28 . 

For the numerical calculation, we list the other input parameters in Table I. 

Using the wave functions and the values of relevant input parameters, we find the 
numerical values of the corresponding form factors — > fo(ss) at zero meomentum 
transfer 

F B^ Mm{q2 = Q) = Q 33+0.024.0.024-0.0^ ^ (2g) 

F f^M™»tf = 0) = -0.2512:8^8^, scenario I, (27) 
F B^M± m{q2 = Q) = .59±8:81±8:gl±8:81, scenario II, (28) 

where the uncertainties are from the decay constant, the Gegenbauer moments Bi and 
B 3 of the meson f Q . The large form factors result from the large decay constants of the 
scalar mesons. The opposite sign of the B° s — > / (1500) (ss) form factor in the upper two 
scenarios arises from the decay constant of /o(1500). These values agree well with those 
as given in Ref. [29 



In the .Barest frame, the decay rate of B® — > foK° can be expressed as 

' G " \M\ 2 (l-r%), (29) 



327rms 

where r/ = rrif /mB s and Ai is the decay amplitude of B — > foK°, which has been given 
in equation (Till) . If / (980) and / (1500) are purely composed of nfi(ss), the branching 
ratios of B° s f (980)K°, f (l500)K° are 



B(B° S ^ f (980)(nn 


K°) 


- 44-0.7+1.2+0.1\ v 
— l°-^-0. 5-1.0-0. OJ X 


KT 6 , 


scenario 


I- 


(30) 


B(B° S ^ f (1500)(nn 


K°) 


- (A A+0.1+3.4+0.5N 
~~ 1^-^-0.1-2.0-0.3/ A 


KT 6 , 


scenario 


I- 


(31) 


B(B° S ^ f (1500)(nn 


K°) 


- (Q Q+0.8+1.4+0.2N 


10~ 6 , 


scenario 


II; 


(32) 


S(S s °^/ (980)« 


K°) 


— l-L-U-O.l-O^-O.lJ X 


10~ 6 , 


scenario 


I: 


(33) 


B(B° S ^ f {1500)(ss] 


K°) 


(n 9 +0.00+0.30+0.06\ 
— l u - z -0.00-0.13-0.05J 


x 10" 


~ 6 , scenario I, 


(34) 


B(B° S / (1500)( S r 


K°) 


_ (n -i +0.5+0.5+0.7^, v 
— ^■ 1 -0.4-0.4-0.5/ A 


10~ 6 , 


scenario 


II. 


(35) 



where the uncertainties are from the same quantities as above. 

In table II, we list values of the factorizable and non-factorizable amplitudes from the 
emission and annihilation topology diagrams of the decays B® — > f (980)K, f Q (1500)K. 
^eL) an d are the K emission (annihilation) factorizable contributions and non- 

factorizable contributions from penguin operators respectively. Similarly, F^ and M^?\ 
denote as the contributions from f emission (annihilation) factorizable contributions 
and non-factorizable contributions from penguin operators respectively. It is easy to see 
that Ffy-> and are larger than F^ and that is fo- emission diagrams give 

large contributions. M/ 0,T denotes the fo emission non-factorizable contribution from tree 
operator 2 . From the table, one can find that the contributions from nn component of fo 
are larger than those from ss component, the one reason is that M/°' T from the operator 
2 is much larger than other amplitudes. Furthermore, B(B® — > fo(nn)K) is about two 
or five times of B(B® -»■ f (ss)K). Certainly, as to the decay B Q S ->■ / (1500)if in SI, the 
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FIG. 2: The dependence of the branching ratios for B° s -> / (980)i^° (a) and B° s -> / (1500)K° 
(b) on the mixing angle 9 using the inputs derived from QCD sum rules. For the right panel, the 
dashed (solid) curve is plotted in scenario I (II). The vertical bands show two possible ranges of 
9: 25° < 9 < 40° and 140° < 9 < 165°. 



TABLE II: Decay amplitudes for 5° -»• f (980)K°, f {1500)K° (xl(T 2 GeV 3 ). 



ss 


F K 

± e 


Mf 




F K 


Mi° 


/ (980)K° (SI) 


-1.0 


0.2 + 0.5* 


-0.4-0.1* 


2.2 - 10.3* 


-5.7-5.1* 


/o(1500)K° (SI) 


1.8 


-0.6 + 0.8*' 


-0.1 - 0.2* 


6.6 + 13.4* 


-2.2 - 8.1*' 


/o(1500)K° (SII) 


-2.4 


0.6 + 0.07* 


-0.5 - 0.04*' 


-0.8 - 21.1* 


-6.3-1.7* 


fin 


Ft 


Mi ' 1 ' 






pfo 

r a 


f (980)(nn)K° (SI) 


7.2 


57.5 + 47.6*' 


-10.1 - 9.5i 


0.4 + 0.4* 


1.9-8.7*' 


/o(1500) (nn)K° (SI) 


-7.6 


20.9 + 77.1*' 


-2.6 - 15.1* 


-0.7 + 0.7*' 


6.8 + H.8* 


f {l500)(nn)K° (SII) 


14.6 


57.9 + 14.3* 


-10.7- 2. 9* 


0.9 + 0.1* 


-1.4-17.8* 



difference is greater. One can see the values of B{B S — > fo (1 500) (nn)K) in SI and SII are 
close to each other, but the values of B(B® — > / (1500) (ss)if) in the two scenarios are 
much differ, which results the total branching ratios have an apparent difference between 
these two scenarios (shown in figur 2(b)). 

In figure [2j we plot the branching ratios as functions of the mixing angle 9. If the mixing 
angle 9 falls into the range of 25° < 9 < 40°, the branching ratio of B® — > f (980)K is: 

2.0 x 10" 6 < B(B° S -> / (980)if) < 2.6 x 10" 6 , (36) 

while 9 lies in the range of 140° < 9 < 165°, B(B° S -> / (980)iT) is roughly 6.5 x 10~ 7 . 
The dependence of the branch ratio of B® — > /o (980) if is strong in the whole mixing angle 
range, but not sensitive in some ranges, for example, 140° < 9 < 165°. As to the decay 
B° — > /o( 1500) if, because there are more discrepancies for the structure of / (1500), we 
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do not show the possible allowed mixing angle range in Fig. 2(b). Lattice QCD predicted 



so 



that the mass of the ground state scalar glueball is around 1.5 — 1.8 GeV |30l . 131 
the three mesons / (1370), /o(1500) and / (1710) become the potential candidates. The 
mixing matrix can be written as 

E3 



/o(1710)\ / Ql a 2 a 3 \ / G \ 

/o(1500) = [h b 2 b 3 \ [ ss . (37) 
/o(1370)/ \ci c 2 cj \nn) 

For each physical scalar meson, the corresponding component coefficients satisfy the nor- 
malization condition, so we have a/ | b\ \ 2 + | b 2 \ 2 + | b^ \ 2 = 1 for the meson / (1500). For the 
earlier lattice calculations predicting the scalar glueball mass to be about 1550 MeV and 
the decay width of /(1500) being not compatible with a simple qq state [33], Amsler and 
Close claimed that / (1500) is primarily a scalar glueball [H]]. But in the SU(3) symmetry 
limit, Cheng et al. reanlyze all existing experimental data and find that /o(1500) is a pure 
SU(3) octet with a ver y ti ny glueball content. Their results for the mixing coefficients are 
given in the following [35J 

7o(1710)\ / 0.93 0.17 0.32 \ / G \ 
/ (1500) = -0.03 0.84 -0.54 ) [ ss . (38) 
./o(1370) / \ -0.36 0.52 0.78 / \nn J 

From above Equation, it is ease to see that /q(1710) is composed primarily of a scalar 
glueball. This conclusion is also supported by an improved LQCD calculation (3l[, which 
predicts that the mass of the scalar glueball is 1710 ±50 ±80 MeV. Here we use the mixing 
coefficients for / (1500) given in Eq. f[3"8"j) and neglect the tiny component of glueball, that 
is | /o(1500)) = 0.84 | ss) — 0.54 | nn). The branching ratios in two scenarios given as: 

B(B° ->■ f Q (1500)K°) = 9.7 x 10~ 7 , Scenario I, 

B(B° S -»■ /o(1500)if ) = 9.5 x 10" 7 , Scenario II, (39) 

which is obtained by the mixing angle taken as —32.7°. From these results, it is difficult 
to determine which scenario is more preferable than the other. 

Now we turn to the evaluations of the CP-violating asymmetries of the considered 
decays in PQCD approach. For the neutral decays — » / (980)i£s, /o(1500)i^5, there 
are both direct CP asymmetry Aq t p and mixing-induced CP asymmetry A™™ . The time 
dependent CP asymmetry of B® decay into a CP eigenstate / is defined as 

A CP (t) = A d S r P (B° s -+ f) cos(Am s ) + A^(B° S -+ f) sin(Am s ), (40) 

with 

J&W -+ f) = ^cP(B° s /) = ^jj, (41) 

\ = r, e -W A ( B ° ^ 11 (42) 
ve A(B® ->• /), {A2) 

where rj = ±1 depends on the CP eigenvalue of /, Am s is the mass difference of the two 
neutral B s meson eigenstates. Here we only give the direct CP asymmetries. 
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TABLE III: Direct CP asymmetries (in units of %) of #° -> f (980)K°, / (1500)if decays for 
nn and ss components, respectively. 



Channel 


Scenario I 


Scenario II 


Bs -> /o(980)(nn)K° 


-57.0 




4°^/ (980)( S s)iC 







5° -> / (1500) (rm)if 


17.6 


-94.8 


5° -> / Q (1500) (ss)if 










FIG. 3: The dependence of the direct CP asymmetries for B® — > /o (980) if (a) and B® 
/o(1500)if (b) on the mixing angle 8. For the right panel, the dashed (solid) curve is plotted in 
scenario I (II). The vertical bands show two possible ranges of 6: 25° < 9 < 40° and 140° < 9 < 
165°. 



In B° s — > fo(ss)K, there is no tree contribution at the leading order, so the direct 
CP asymmetry is naturally zero. As B Q S — > f (nn)K, the corresponding direct CP asym- 
metries are listed in table III. For the decay B° s — > /o(1500)if, the amplitudes from the 
non-factorizable /(remission and annihilation topologies, that is M/° and _F/°, are con- 
structive in scenario II but are destructive in scenario I (seen in table II). Furthermore, the 
contributions from the factorizable /(remission diagrams have an opposite sign between 
scenario I and II (which is because that the decay constant of / (1500) is opposite in 
two scenarios). These reasons result that there exists a great difference for the direct CP 
asymmetries of B° s — > /o(1500)(nn)if in two scenarios (shown in table III). The depen- 
dence of the direct CP violating asymmetries for these decays are shown in Fig.3. From 
Fig. 3 (a), we can see that the signs of the direct CP asymmetries in the two allowed mixing 
angle ranges are opposite, it gives the hint that one can determine the mixing angle by 
comparing with the future experimental results. For the decay B° s — > / (1500)if, if we 
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still use the mixing scheme given by Eq.( l38l for /o(1500), the direct CP asymmetry is 

Acp(B° -»■ f (1500)K°) = 47.4%, scenario I, 

Acp(B° 8 -> / (1500)K°) = 62.8%, scenario II. (43) 

Although the CP asymmetry for the decay B® — > f (1500)K° is large, it is difficult to 
measure it, since its branching ratio is small. 

V. CONCLUSION 

In this paper, we calculate the branching ratios and the direct CP-violating asymme- 
tries of B° s — > /o (980) K, f (1500) K decays in the PQCD factorization approach. From 
our calculations and phenomenological analysis, we find the following results: 

• In general, the contributions from /o-emission diagrams are larger than those 
from if-emission diagrams. Especially, the /(remission non-factorizable contribu- 
tion from tree operator O2 is quite larger than other amplitudes. For the decays 
Bg — » fo(980)K,fo(1500)K, the contributions from nn component are larger than 
those from ss component in two scenarios. 

• Using the wave functions and the values of relevant input parameters, we find the 
numerical results of the corresponding form factors B® — > fo(ss) at zero meomentum 
transfer 

F 5 °-^ 98 V = 0) = 0.331°;^?^, scenario I, (44) 
if° -> /0(150( V = 0) = -0.251JS185K8S, scenario I, (45) 
^ /0(150 °V = 0) = 0.59±S5g±8^±gS, scenario II. (46) 

The values of Fq°^ M1500) (q 2 = 0) for two scenarios can be used to identify which 
scenario is favored by compare with the future experimental results. 

• If the mixing angle 9 falls into the range of 25° < 9 < 40°, the branching ratio of 
B° s ->• /o (980) if is 

2.0 x lO" 6 < B(B° S ->• /o(980)iT) < 2.6 x 10~ 6 , (47) 

while 9 lies in the range of 140° < 9 < 165°, B(B° S ->• / (980)iT) is about 6.5 x lO" 7 . 

• if we identify the meson / (1500) as a pure SU(3) octet state and use the mixing 
scheme giving by | /o(1500)) = 0.84 | ss) — 0.54 | nn), one can find that the 
branching ratios in two scenarios are both close to 1.0 x 10~ 6 . Although the CP 
asymmetry for the decay B° s — > f ( 1500) K° is large, it is difficult to measure it, 
since its branching ratio is small. 
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